Coulombian Disorder in Periodic Systems 
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We study the effect of unscreened charged impurities on periodic systems. We show that the 
long wavelength component of the disorder becomes long ranged and dominates static correlation 
functions. On the other hand, because of the statistical tilt symmetry, dynamical properties such as 
pinning remain unaffected. As a concrete example, we focus on the effect of Coulombian disorder 
generated by charged impurities, on 3D charge density waves with non local elasticity. We calculate 
the x-ray intensity and find that it is identical to the one produced by thermal fluctuations in a 
disorder- free smectic-A. We discuss the consequences of these results for experiments. 

PACS numbers: 



I. INTRODUCTION 

The effect of quenched disorder on various condensed 
elastic systems is one of the fascinating problems in 
statistical mechanics. Examples of physical systems 
range from domain walls in magnetic and ferroelectric 
materials contact lines of of a liquid meniscus on 
a rough substrata^, crack propagation*, to vortex lat- 
tices in type II superconductors^, charge density waves 
(CDW)^'- and Wigner crystala^iiSiiiii^. In these sys- 
tems, the competition between elastic interactions which 
tend to impose some long range order in the system 
and quenched disorder, leads to the formation of glassy 
phases. Two broad classes of elastic systems can be dis- 
tinguished: random manifold systems such as domain 
walls, contact lines and cracks, and periodic systems such 
as vortex lattices, charge density waves, and Wigner crys- 
tals. The latter are characterized by a long range crys- 
talline order in the absence of disorder and thermal fluc- 
tuations. For these systems, a crucial question is whether 
a weak disorder destroys entirely the crystalline order, or 
whether some remnants of the underlying periodic struc- 
ture remain observable. 

One of the earliest attempts to answer this question, 
was the pioneering work by Larkin^'^ on vortex lattices. 
Using a random-force model, he showed that due to the 
relevance of disorder in the renormalization group sense, 
long range order was entirely destroyed below four di- 
mensions. Above four dimensions, long range order per- 
sists as disorder becomes irrelevant. A similar conclu- 
sion was reached by Sham and Patton for the case of 
a CDW with short range elasticitjJ^, where, using an 
Imry-Ma approach, they concluded that long range or- 
der was impossible in the presence of disorder below four 
dimensions. The problem of short range disorder in pe- 
riodic s ystems with short ranged elasticity was revisited 
m Refs. Ilfill7ll8ll9l It was argued that the periodicity 
present in systems like CDW and vortex lattices plays a 



pivotal role in determining the physics of the system in 
the presence of disorder. More precisely, it was shown 
that, though the disorder is relevant below four dimen- 
sions, due to the underlying periodicity of the system a 
quasi long-range order persisted for dimensions between 
two and four. This is in stark contrast to the earlier re- 
sults which predicted a total destruction of order. The 
resulting phase, nicknamed Bragg glass phase, possesses 
both quasi long range order and metastability and glassy 
properties. ^^•"'^^ It was further shown that the Bragg glass 
phase is stable to the formation of defectsii^iSSiSiiSS Re- 
cent neutron scattering experiments on vortex lattices 
have furnished clear evidence for the existence of such a 
phasei^ 

A complication arises in charged periodic systems due 
to the Coulomb repulsion, which renders the elastic- 
ity non-localiS^iSSiSSiSl This non local elasticity tends to 
rigidify the system, so that short range correlated disor- 
der could be irrelevant in dimension smaller than four.— 
For instance, within the random force model, the corre- 
lation function of the displacement in three dimensions 
displays a logarithmic growth indicating quasi-long range 
order. ^^•^'^■^^ In fact, when the periodic structure of the 
CDW is properly taken into account, the growth of the 
displacement correlation function is even weaker, increas- 
ing only as log(log(r)) with the distance rM A second 
complication arising from Coulomb interaction is that the 
disorder induced by charged impurities has long-range 
correlations. This type of disorder can exist in certain 
doped CDW material£^° such as KMoi^xV^Oa. 

In this paper, we study the effect of the competition of 
the non-local elasticity produced by the Coulomb interac- 
tion with the long-range random potential resulting from 
the presence of charged impurities on the statics. The 
paper is organized as follows: in Sec. m we introduce a 
decomposition of the Coulomb potential on the Fourier 
modes of the periodic structure. With this decomposi- 
tion, we show that only the long-wavelength component 
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of the random potential i.e., forward scattering disorder, 
possesses long range correlations. Using statistical tilt 
symmetrji^, we deduce that due to the short ranged 
nature of the backward scattering terms engendered by 
the disorder, the dynamical properties in the presence 
of charged impurities are not qualitatively different from 
those in the presence of neutral short ranged impurities. 
In Sec. mil we consider the problem of the CDW sys- 
tem, and we derive the non-local elastic Hamiltonian. 
In Sec. II VI we derive the static displacement correlation 
functions and x-ray intensity of the CDW with charged 
impurities and we highlight the similitude of the latter to 
the x-ray intensity of smectic-A liquid crystals subjected 
to thermal fluctuationsi^ In Sec. El we discuss the ex- 
perimental significance of our result and suggest that the 
smectic-like correlations should be observable in experi- 
ments on KMoi-xVxOa. Finally, we summarize the pos- 
sible behavior of the static correlators in a pinned Charge 
Density wave according to the local or non-local charac- 
ter of elasticity and the presence or absence of charged 
impurities. 



II. ELASTICITY AND DISORDER IN 
PERIODIC SYSTEMS 

In this section, we discuss how Coulomb interactions 
affect elasticity and disorder in periodic systems. For a 
periodic elastic structure, the density can be written as: 



p(r) = po{r) 



G 



'U(r)) 



(1) 



where po{r) = pq{1 — W ■ u) describes the density fluctua- 
tion arising from the long wavelength deformation of the 
periodic structure and po is the average density. In the 
second term, the vectors G belong to the reciprocal lat- 
tice of the perfect periodic structure, and u(r) represents 
a slowly varying^^ elastic deformation of the structure.-^ 
The quantities e**^'^''""'^'"-'-' describe fluctuations of the 
density on the scale of a lattice spacing. The low en- 
ergy physics of the periodic structure can be described 
in terms of a purely elastic Hamiltonian which has the 
generic form for isotropic systems 



(2) 



where c is the elastic coefficient and is a shorthand for 
J (Pr. This form can easily be generalized to anisotropic 
systems. Well known examples of charged periodic struc- 
tures are the Wigner crystal)2ii2iiiiSi24 charged colloidal 
crystalsji^^ and the charge density wavesi^S'^LSS In many 
charged systems, unscreened Coulomb interactions are 
present: 



P(r)p(r') 
47re|r — r'l 



(3) 



and strongly affect the elasticity and dispersion of the 
compression modes of the system. Moreover, in the pres- 
ence of charged impurities, the original charge density on 
the lattice interacts with the charge impurity yielding: 



Hdis. — e 



P(r)pimp.(r0 
47re|r — r'l 



(4) 



where pimp, denote the impurity density. Using the 
decomposition of the density 1^, we now show that 
the Coulomb interactions fundamentally modify only the 
long wavelength components of the elasticity and of the 
disorder energy. 

To better handle the periodicity of the elastic struc- 
ture, it is convenient to use the decomposition of the 
Coulomb interaction in terms of the reciprocal lattice 
vectors G. In three dimensions, this decomposition reads 



1 



47r|r| J (27r)3 q2 



G 



where 



^G(r) 



d3 



q e 



Bz (27r)3 (q + G)2' 



(5) 



(6) 



and Jg^ indicates that the integral is restricted to the 
first Brillouin zone. It is straightforward to check that 
V-dr) — V^{r). Using Eq. (0), the interaction term He 
can be rewritten as: 



He ^ -J2 f Vb(r-r')e'°-("("^)-"("-'» 

+^^^/o(r-r>o(r)po(r'). (7) 



Note that due to the slow variation of u(r), terms involv- 
ing the oscillatory factors e*(*^~*^ ^''^ can be dropped from 
the interaction. Let us first consider the term involving 
long wavelength fluctuations of the density. Since we are 
interested only in the long wavelength properties, we can 
replace the integration over the Brillouin zone in Vo(r). 
by a Gaussian integration: 



f f 

Jbz (2^)3 ^ J 



{2n) 



(8) 



with the parameter a chosen so that tt/o \Gy- 



being the reciprocal lattice vector having the 



shortest length. In this case, Vb(r) can be obtained indi- 
rectly by solving the Poisson equation with a Gaussian 
charge density and is found to be 



Voir) 



1 



-erf 



(2a) ' 



(9) 



In the limit r ^ a, we recover the known result Vq ~ 
l/(47rr). Clearly, the non-oscillating component of the 
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Coulomb potential remains long-ranged and tends to 
rigidify the system. 

It now remains to be seen whether the oscillating parts 
of the Coulomb interaction specified by Vq for G are 
long-ranged or not. We first note that the above trick 
of replacing the integration over the Brillouin zone by a 
Gaussian integral over the entire space is not applicable 
anymore, as it would introduce a spurious integration 
over a region where G -|- q = 0. This would result in an 
(incorrect) 1 jr behavior of Vq^o ('') ■ To obtain a correct 
estimate for Vq we replace the integral over the Brillouin 
zone by 

f d^q f d^q , , 

where FBz{q, e) is an indefinitely derivable hmction with 
a compact support contained in the first Brillouin zone 
(see App. ^ for an explicit form of -Fbz(<Z, e))-^'^ Obvi- 
ously, for G 7^ 0, FBz{q, e)/|q+Gp is also an indefinitely 
differentiable function of compact support. A well-known 
theorem2M£ then shows that the Fourier transform of 
FBz{q,£)/\q+ is indefinitely differentiable and for 
r ^ oo is o(l/r"') for any n > 0. This implies that the 
function Voir) is short ranged. Incorporating the above 
results in Eq.(|Zj), we see that while the non oscillating 
part of the Coulomb interaction modifies the long wave- 
length behavior of the elasticity, rendering it non-local, 
the short ranged nature of the oscillatory terms merely 
renormalizes the elastic coefficients. This is explicitly 
shown in Add. iBlfor the particular case of a CDW. The 
resulting non-local character of the elastic interactions 
modifies strongly the static and dynamic properties of 
the system,SLaii4^ 

To understand the nature of the interaction with the 
charged impurities, we use the above procedure to rewrite 
the random potential generated by the impurities as: 

, . ^ f d\ e'q--p,„p.(q) 

= ^e^«-f/G(r). (11) 

G 

Using this in Eq.Q, the interaction of the system with 
the random potential is given by 

^ f t/G(r)e»«"("-) 

+- / ;7o(r)po(r), (12) 

In Eq. (|12|l . the interaction of po with the ran- 
dom potential Uo is called forward scattering, and 
the terms containing e'*^'"'^''^ are called backward scat- 
tering. This nomenclature originates in the the- 
ory of electrons in ID random potentiali^ To cal- 
culate the disorder correlation functions, we consider 
the case of Gaussian distributed impurities where^"* 



Pimp.(G + q)pi„,p.(G' + q') = (27r)3D(5G,-G'5(q + q'), 
the parameter D measuring the disorder strength. Con- 
sequently, we find that for G 7^ 0: 

Using the same arguments as before, we infer that the 
correlations of f/G(r) are short ranged, as in the case of 
neutral impurities, for all G except G = 0. This im- 
plies that the backward scattering terms induced by dis- 
order are short-ranged and the treatment of these terms 
within the replica or the Martin-Siggia-RoseiSii& meth- 
ods is identical to the case of neutral or screened impu- 
rities. However, the G = component 

manifests power law decay of the forward scattering cor- 
relations. This term however can be gauged out by the 
statistical tilt symmetry^, and affects mainly the static 
properties of the periodic system. Typically, in periodic 
systems with both short range disorder and local elas- 
ticity, the contribution of the forward scattering disorder 
can be neglected and it is the backward scattering that 
induces collective pinning and Bragg glass features like 
a quasi order in the static correlation functions. Here, 
we have shown that even in the case of long range dis- 
order, the backward scattering terms behave essentially 
like their short ranged (neutral impurities) counterparts. 
However, the effect of the forward scattering terms on 
the correlation functions has to be studied carefully. In 
the next section, we show that in the case of charged 
impurities in a charge density wave system, the forward 
scattering term strongly modifies the static correlation. 
Finally, we remark that our decomposition of the elastic 
energy and the impurity potential is not exclusive to the 
Coulomb potential and is applicable to other long-range 
potentials. As a result, the conclusions of the present 
sections are expected to be valid for more general long 
range potentials. 

III. CHARGE DENSITY WAVES 

In this section, we re-derive the elastic Hamiltonian 
for a d dimensional CDW with screened Coulomb inter- 
actions between the density fiuctuations at zero temper- 
ature. We consider an incommensurate CDW, in which 
the electron density is modulated by a modulation vector 
Q incommensurate with the underlying crystal lattice. In 
this phase, the electron density has the following form^: 

P(r) = Po + • V(A(r) + pi cos(Q • r + <^(r))) , (15) 

where po, is the average electronic density (see App. El 
for details). The second term in Eq. H15|) is the long 
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wavelength density and corresponds to variations of the 

density over scales larger than . The last oscillating 

term describes the sinusoidal deformation of the density 

at a scale of the order of Q"^ induced by the formation 

of the CDW with amplitude pi and phase (p. 

In the absence of Coulomb interactions, the low energy 

properties of the CDW can be described by an effective 

Hamiltonian for phase fluctuations. For CDW aligned 

along the x axis, i.e., Q = Qx, this phase only Haniilto- 
nian reada47,48,49,50. 



An 



V 



(16) 



where is the Fermi velocity and ric is the number 
of chains per unit surface that crosses a plane orthog- 
onal to Q. The velocity of the phason excitations par- 
allel to Q is Vx — {irLe/m^y/^VF with m* the effective 
mass of the CDW and the mass of an electron. Vy 
and Vz denote the phason velocities in the transverse di- 
rections. A crucial observation is that a deformation of 
(j){r) along Q produces an imbalance of the electronic 
charge density which then augments the electrostatic en- 
ergy due to Coulomb repulsion between density fluctua- 
tions. We evaluate the contribution of Coulomb interac- 
tions screened beyond the characteristic length A which 
accounts for the presence of free carriers. This length 
diverges in the limit T O.^^'^^ The electrostatic energy 
takes the form: 



Hr 



e 

8^ 



d'^vd'^v'e 



(17) 



where we have assumed for simplicity an isotropic dielec- 
tric permittivity e of the host mediumi^iiSiSiiSi^ Due 
to the periodicity of the CDW system, we can use the de- 
composition of the Coulomb potential derived in SecHll 
obtaining: 



He 



2e 



r,r' 



(18) 



r - r' 



VLQ(r-r')e*('^('')"'^('"'» + 



c.c. 



In Eq. H18|l . we have neglected the contribution of the 
higher harmonics of the CDW. Note that the oscillating 
terms, as discussed in App.0 only contribute to a renor- 
malization of the coefficients in the short range elastic 
Hamiltonian H16|l and thus can be neglected. However, 
the contribution of the long-wavelength term has more 
dramatic effects and reads 



He 



2eQ2 



BZ 



A-2 



(19) 



in the three dimensional case. It is interesting to note 
that Coulomb interactions generate a non local elasticity 
i.e., a q-dispersion in the elastic constant. 



The total Hamiltonian in d = 3 now reads. 



H 



cl. 



(20) 



G-Hq) - 



2n 



y ^2 



2 „2 



(q2 + x-^)e 



1x + -^1y + —1. 



where the lengthscale ^ is defined by: 



Hchvp 



27re2p^ 



(21) 



Depending on the ratio A/^, two regimes of behavior can 
be identified: (i) Short ranged elasticity: when A/^ ^ 1 
the Coulomb correction to the short range elasticity is 
small even in the limit g — )■ and hence can be neglected, 
(ii) Long range elasticity: for A/^ ^ 1, the Coulom- 
bian correction to the short-range elasticity cannot be 
neglected. This regime is relevant at low temperatures, 
when the number of free carriers available to screen the 
Coulomb interaction is suppressed by the CDW gapMi^ 
Mean field calculations show that this regime is obtained 
for temperatures T < 0.2Tc where Tc is the Peierls transi- 
tion temperature^ In the following, we focus on regime 
(ii), and accordingly, we take A~^ = in Eq. H2()(l . 



IV. FORWARD SCATTERING 

As discussed in SecHJ the case of the short ranged elas- 
ticity has been studied by various authors. For charged 
periodic systems with short range disorder and a non- 
local elasticity generated by Coulomb interactions, it is 
known that d — S becomes the upper critical dimen- 
sion for disorder and the displacement correlations grow 
as B{r) = loglog AriSSiSS Here, we study the effect of 
the long range disorder on the static correlations of a 
charged periodic system. Since, the backward scattering 
terms generated by such a disorder are short ranged, they 
lead to the same physics as that of short ranged disorder 
with the corresponding nonlocal elasticity. These terms 
contribute a log log r term to the displacement correla- 
tions. However, in this case a simple dimensional analy- 
sis shows that the forward scattering terms generate the 
leading contribution to the correlation functions. In the 
following, we calculate the contribution of the forward 
scattering disorder to the displacement correlation func- 
tion in the d = 3 CDW. 



A. Displacement correlation functions 

The displacement correlation function is defined by: 



B{r) = mr)~mr) 

= ;^ E WM^n - cosq • r] . (22) 
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The calculation of the correlation induced by the for- 
ward scattering disorder is analogous to the calculation 
of Larkin for the random force modelii^ Assuming an 
infinite screening length A, the Hamiltonian reads: 



where ^' = ^v± /vx, r\ = + and is a momentum 
cut-off. A study of the limits of this function for — > oo 
and jxl ^ cx) shows that its asymptotic behavior is well 
described by: 



el. 



Aire 



(23) 



Using Eq. H15() in Eq. (|23|l . we obtain an expression of 
the form Eq. H12|) . Keeping only the forward scattering 
term we get: 



H = 



(Pq 



—Pirap.{-q)(P\<l) 



Qeq^ 



Shifting the field (j) 



e^po iqxGjq) 
Qe g2 



(25) 



brings the Hamiltonian (|24|l back to the form of Eq. (|20|l . 

The average over disorder now yields: 



+ 



Pimp, (q) Pimp, {-q) 



TG{q) 



e'pl qlGjqy 
(52 £2 qi 



D 



(26) 



where (...) and denote thermal average and disor- 
der average respectively. Eq shows that even in 
the presence of Coulombian disorder, the statistical tilt 
symmetry^ is preserved. This implies that in the pres- 
ence of backward scattering disorder, the forward scatter- 
ing term can be gauged out by H25|) , and the contribution 
of the forward scattering disorder B^^ is simply added to 
the one obtained from the backward scattering disorder, 
5BS,22^^g conclude 



B^^{r) = 2D 



d\ qlG{qf 



Q^e^J (2^) 



q^ 



[l-cos(q.r)]. (27) 



We want to evaluate this integral for the case of Vy = 
Vz — v±. In the following we will use qj^ — qy + q'^. To 
obtain the asymptotic behavior of B(r) for r — > 00 we 
need to consider the q ^ limit of the integrand. The 
form of G{q) suggests a scaling q^ ~ which then allows 
us to consider the integral: 



Fir) 



dh 



qx 



(1 - cos(q • r)) (28) 

„2 



(27r)3 [ql + i^'qlf]^ 



F{v) 



■In 



A7^ 



Therefore, we have for r — > 00: 

'ri +4(w^|x|CM) 



B^^(r) = Kin 



A^ 



where: 



(24) 



DQ^Vx 



(29) 



(30) 



(31) 



The full asymptotic correlation function is given by the 
sum of the forward scattering contribution, Eq. (|3()|l . and 
the backward scattering contribution given in Eq. (51) of 
Ref. J£ for the case of a short-range disorder and non- 
local elasticity: 

B^^ir) = log [log imax{A\x\, (Ar^)^))). (32) 

Obviously, the contribution of the backward scattering 
terms is subdominant and can be neglected. 



B. Analogy with smectics-A 

We note that the result Eq. 1)29(1 can be obtained in 

the entirely different context of liquid crystals. If we 

consider a smectic-A liquid crystal, its elastic free energy 
reada57,58,59,60. 



(33) 



where u represents the displacement of the smectic layers, 
B is the compressibility, and fcn measures the bending 
energy of the smectic layers. If we now assume a random 
compression force given by: 



d^rr]{r)dzu{T), 



ri{r)ri{r') = DS{r — r'). 



(34) 
(35) 



a straightforward calculation shows that the displace- 
ment correlation function (u(r) — u(0))2 is given by 
Eq. (|29|l . Smectics-A with disorder have been consid- 
ered previously in Ref. |^ albeit with a different type of 
disorder couphng to \7±u. This yields a displacement 
correlation function superficially similar to F{r) with q\ 
replacing q^ in the numerator. The random compression 
force, which is not natural in the smectic-A context, is 
thus easily realized with charge density wave systems. 
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V. EXPERIMENTAL IMPLICATIONS 

In the preceding sections, we have shown that the for- 
ward scattering terms generated by charged impurities 
lead to smectic-hke order in a charge density wave ma- 
terial. A frequently used technique to characterize posi- 
tional correlations in CDW systems is x-ray diffraction.^^ 
In the present section, we provide a calculation of the x- 
ray intensity resulting from such a smectic-like order, and 
we provide a quantitative estimate of the exponent k. 



A. x-ray intensity 

The intensity of the x-ray spectrum is given byS^ 

= T^T. e-^'(^'-^^-) (/J,e-»?(«--".)) . (36) 



Ui is the atom displacement from the equilibrium position 
Ri, fi represents the total amplitude scattered by the 
atom at the position i and depends exclusively on the 
atom type. We consider the simple case of a disordered 
crystal, made of one kind of atoms, characterized by the 
scattering factor / — A//2, and containing impurities of 
scattering factor / -I- A//2. Since we are interested in the 
behavior of the scattering intensity near a Bragg peak 
{q ^ K), we can use the continuum approximation.^^ In 
the case of the CDW, the lattice modulation is given by: 



u{t) = [cos((3a; + (?!)(r))] 



(37) 



It is well known that the presence of a CDW in the com- 
pound is associated with the appearance of two asymmet- 
ric satellites at positions q ^ K zL Q around each Bragg 
peakiS The intensity profiles of these satellites give ac- 
cess to the structural properties of the CDW. For this rea- 
son a lot of work has been done to compute and measure 
these intensitiesi22i2Si^SiS4iS^ By expanding Eq. for 
low q{ui — Uj), one finds an expression of the x-ray satel- 
lite intensity comprising a part Id, which is symmetric 
under inversion around the Bragg vector K and a part Ja 
which is antisymmetric under the same transformation^ 
The symmetric part is given by the following correlation 
function: 



/d(q) = /Vy d 



u{r/2)ui~r/2)), (38) 



and the antisymmetric part by: 

4(q) 
a • (b X c) 



- 2gA/Im ^ d^r e-'^i-(pi„p(r/2)^x(-r/2)), 



where: Sq ^ {q — K) ~ Q, and a • (b x c) is the volume 
of the unit cell of the crystal. After some manipulations. 



Eq. H38|) can be rewritten as: 

h{K + Q + k) = ulJ^K^ J dhe-^'^ '^Cdir), (40) 

UK + Q + k) = -jKuoAfVJVD J d^r e-^'' ''Ca(r) 

where Af is the number of impurities in the unit cell, and: 
Cd(r) <Je'«'-/2)-0(-r/2))^ ^ (41) 



a(r) = x(r)Cd(r), 



(42) 
(43) 



where x(r) is defined by Eq. (33) of Ref. 29,. It is easy 
to show, using this definition and the statistical tilt sym- 
metry that x(r) is independent of the forward scattering 
disorder. In Eq. H41|l . Cf '^ is the backward scattering 
contribution which has been obtained in Ref. |2^ and 
Cj'^ is the forward scattering contribution, given by: 



A- 



(44) 



where we have used Eq. (|30(l . assuming a Gaussian dis- 
order. Using Eq. (|32(l . one sees that the term C^'^ gives 
only a logarithmic correction to H41|) . As a result, the 
symmetric structure factor 1^ is dominated by the con- 
tribution of the forward scattering disorder. To obtain 
the structure factor, we Fourier transform Eq. (|44|l to 
obtain 



+ 00 



JO {l + u)-f 
Using the following relation. 



du f\qx\r'ivx 
cos — — u 



+ 00 



(l + u) 



■ cos(Au) 



r(7) Jo 



— 7(45) 



we finally obtain 



/d(q) = 



22(,.-i)r(K) {£,vl/vxY-^ 



dw 



w 



1 



(46) 



so that /d(q) ^ {\qx\Y ^ for lliiv^/vx) < |fe| and 
^d(q) ^ (Igxl)^*"^^^ otherwise. The intensity /d(q ^ 0) 
is divergent for for k < 2 but is finite for k > 2, i.e. for 
strong disorder. Next, we turn to the evaluation of /a. 
From Ref. |2^ we know that xl'') ^ when a;^ ^ r\ 
(39) and x('') ~ l/''"! when ^ r\. This implies that 



la is subdominant in comparison with I^i. In particular, 
4(q) {\<lx\T~^ for q\_{ivi^/vx) < \qx\ and /a(q) ~ 
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(|(7_L |)^("~^^ otherwise. We illustrate the behavior of the 
x-ray intensities on Fig. ^ 

We note that these intensities are remarkably similar 
to those of a disorder-free smectic-A liquid crystatfS, at 
positive temperature. In fact, the expression of the ex- 
ponent K Eq. H31|l is analogous to the expression (5.3.12) 
m Ref . fsgl, with the disorder strength D playing the role 
of the temperature ksT m the smectic-A liquid crystal. 





K-Q 



K+Q 



FIG. 1: A sketch of the x-ray intensity in a CDW with 
Coulomb elasticity and charged impurities for K, Q parallel to 
the chain direction. We have taken k = 0.5 in the expressions 
of la and Id- 



Elasticity 


Short range 


Disorder 
Long 


range 


Local 
Non-local 


Unphysical 


Hi) 
Id{qx) 
I4<1±) - 


|g^|2(«-2) 



TABLE L The different cases with short range and long range 
random potential and elasticity. 77 ~ 1 is the Bragg Glass 
exponenli— . k is defined in Eq. (1311 . 



Kq.sMoi-xVxOs, we find that the disorder strength can 
be expressed as a function of the doping and obtain: 



D = x{l 



#(Mo atoms/unit cell) 
a • (b X c) 



(51) 



This formula is derived in App. El For the crystal pa- 
rameters, a =18.25A, b =7.56A, c =9.86A, /3 =117.53°^, 
with 20 Molybdenum atoms per unit cell, and a dop- 
ing X = 3%, the disorder strength £> = 4.8 x lO^^m""^. 
Moreover, using the experimental bounds of the veloc- 
ities: 3.6 X lO^m.s^^ < t;^ < 1.6 x lO'^m.s"-^ and 
Vx ~ 3.7 X lO'^m.s"^, we find that k is in the range 
[0.16 — 0.8]. Therefore, the smcctic-like order should be 
observable in x-ray diffraction measurements on this ma- 
terial. 



B. Estimate of the exponent k. 

Let us turn to an estimate of the exponent k appearing 
in the intensities to determine whether such smectic-like 
intensities are indeed observable in experiments. To do 
this, we first need to determine whether Coulomb inter- 
actions are unscreened by comparing the screening length 
with ^ given by Eq. (|21|l . This question is relevant only 
to a material with a full gap, in which free uncondensed 
electrons cannot screen charged impurities. A good can- 
didate is the blue bronze material K0.3M0O3 which has 
a full gap, and is well characterized experimentally. We 
now evaluate the quan tity ^ for this material. Using the 
parameters of Ref. |6y: 

He = 102°chains/m^ (47) 

VF = 1.3 X 10V.s~\ (48) 

po = 3 X lO^^e^/m^^ (49) 

g = 6 X l0^m-\ (50) 

and a relative permittivity of cko 3M0O3 = Ij so that e 
in Eq. H21|l is equal to the permittivity of the vacuum, 
we obtain: ^ ~ 5A. Therefore, the screening length can 
be large compared to ^ at low temperature, and we ex- 
pect that Coulombian effects will play an important role 
in this material. We can use this value of ^ to evalu- 
ate the exponent k in Eq. H31|) . For the doped material 



VI. CONCLUSION 

In this paper, we have introduced a decomposition of 
the disorder induced by charged impurities in terms of 
the reciprocal lattice vectors of a periodic charged elastic 
system. Using this decomposition, we have shown that 
only the long wavelength (forward scattering) component 
of the disorder was long-range correlated. Components 
with wavevectors commensurate with the reciprocal lat- 
tice of the periodic elastic system remain short ranged. 
The latter can thus be treated with the standard tech- 
niques developed for impurities producing short range 
forcesii^ We find that only the forward scattering is af- 
fected by the long-range character of the forces created 
by charged impurities. Due to the statistical tilt sym- 
metry, this implies that only the statics of the periodic 
elastic system is modified by Coulombian disorder. This 
has allowed us to obtain a full picture of the statics of 
charge density wave systems in d = 3 in the presence of 
charged and neutral impurities. The results are summa- 
rized in tabled A remarkable result is that in the case of 
charged impurities in a system with unscreened Coulomb 
elasticity, the x-ray intensity turns out to be identical to 
that produced by thermal fluctuations in a smectic-A liq- 
uid crystals, with the disorder strength playing the role 
of an effective temperature. This behavior of the scat- 
tering intensity should be observable in the blue bronze 
material K0.3M0O3 doped with charged impurities such 
as Vanadium. 
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APPENDIX A: DECOMPOSITION OF THE 
COULOMB POTENTIAL 

In this appendix, we give more details on the decom- 
position of the Coulomb potential using infinitely difFer- 
entiable functions of compact support. First, let us dis- 
cuss a simple decomposition in ID. We consider the func- 
tion F^{x) such that: 



F,{x) = F,{-x) 
F^{x) = 1, < X > 1 

If , r 2(1 --x) 

F,{x) = 0,x>l + e 



l\<e 
(Al) 



It is easy to check that is continuous, infinitely differ- 
entiable, and that: 



E 



F,{x - 2n) = 1. 



(A2) 



Applying this formula to a one dimensional reciprocal 
space, we obtain: 



nx — — oo 



(A3) 




-2n/a 



-n/a 



n/a 



2%la 



FIG. 2; The graph of the function Fe^qa/ir), of compact sup- 
port for a value of e = ... (solid line). Dotted lines repre- 
sent Fe{qa/n ± 2). We can graphically check the validity of 
Eq. lIXsl . 

2n/Q 




i.e. we have constructed explicitly a partition of the 
unity.^^ The generalization to a cubic lattice in a three 
dimensional space is obvious: 



oo 



nx,ny,nz— — oc 



X 



FA- 



ly 



xF, ( - 

TT 



27rri, 



= J2Frici-G) = l. 



G 



FIG. 3: The ground state of an array of CDWs coupled by a 
repulsive interaction. 



wave with unscreened Coulomb interactions, Eq. H17|l . 
and show that they only induce corrections to the short 
range elasticity. Inserting the expression of the density, 
Eq. ((T3|l . in Eq. Q, the contribution of the oscillating 
component of wavevector Q is given by, 



(A4) 



The function F^^ has a compact support, and vanishes 
rapidly outside the Brillouin zone. 



rjosc. 



9 9 

e P i 
2e ^ 

n,n' 



^ J dxdx'VQ{x — a;', n — n') 
X cos[(?!)(x, n) — (/)(a::', n')], (Bl) 



APPENDIX B: CONTRIBUTION OF THE 
OSCILLATING COMPONENTS OF THE 
DENSITY TO THE ELASTIC HAMILTONIAN IN 
THE PRESENCE OF COULOMB INTERACTION 

In this appendix, we calculate the contribution of the 
oscillating terms to the Hamiltonian of the charge density 



where we have reestablished the discrete character of the 
transverse dimension y. Both the intrachain (n ~ n') 
and the interchain contributions (n ^ n') arc short 
ranged. Let us first consider the case of n ^ n'. We 
have to compute integrals: 



Vqix.n) 



gi(gxa;+qi-n) 



(B2) 
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To evaluate the above integral in closed form, we need 
to make some approximations. Since vanishes for 
Q + Ix — 0, we can neglect compared to Q in this 
integral. Then, we can extend the integration over the 
whole reciprocal space without encountering any singu- 
larity. The Qx integration produces a S{x) function, and 
the q± integration gives an exponential in 2D and a mod- 
ified Bessel function Kq in 3D. In the two-dimensional 
case, we find that: 



where po is the average electron density, pi is the con- 
densate amplitude at T = 0, takes into account the 
reduction of CDW order by thermal fluctuations {\ip\ — 1 
at T = 0), the factor jjc takes into account the pres- 
ence of non-condensed electrons^^ at finite temperature 
(at T ^ 0, pc ~ 1) and Q ~ 2kp. Using the relation 
kp = fPO: valid in a one-dimensional system, one can 
see that this relation simplifies (at T = 0) to: 



VQix,n) 



2Q 



and in the three-dimensional case: 



Ko (Q VM)^ + MV) 



2tt 



(B3) 



5{x). (B4) 



where, iy and are interchain spacings. Due to the 
exponential decay of the interchain interaction with the 
distance, it is justified to neglect interchain interactions 
beyond nearest neighbors. The logarithmic divergence in 
Eq. (|B4p for n = is an artefact of the approximation we 
make when we integrate over the entire reciprocal space 
instead of the first Brillouin zone. A more refined esti- 
mate yields a finite, short ranged intrachain contribution. 
The short range contribution in the electrostatic energy 
thus reads: 



TTOSC. 



E 



dx cos[4){x, n + Ba) — (j}{x, n)] 



J Vq{x — x' ,0)dxdx' cos[(j){x,n) 



where: 



2 2 

Zne 



(/)(a:,n)], 
(B5) 

(B6) 



In Eq. (|B5p . we expand cos(0(a;,n) — (j){x',n)) ^ 1 — 
{x — x')'^ /2{dx4>{x,n))'^ to show that the backscatter- 
ing term reduces to short ranged elastic forces. Making 
(j){x,n) — 0(2:, n) -I- HyTT -I- UzTT, we can make the sign 
in front of Jq negative, and obtain the ground state for 
^ = 0. In its ground state, the CDW is out of phase on 
two nearest neighbor chains. This ground state is repre- 
sented on Fig. 121 The excitations above this ground state 
are described by the Lagrangian ljD3|) derived in App.IdI 



APPENDIX C: ZERO TEMPERATURE LIMIT OF 
THE CHARGE DENSITY IN A CDW 

In the present appendix, we discuss the zero temper- 
ature limit of the expression of the charge density in a 
CDW. Let us consider the form of the charge density in 
the presence of a non-uniform (f>. It is givenSSi^LZSiZi by 
the expression: 

p(r) = PO + ^9x<l>{r) + pM cos{Qx + 0(r)) (CI) 



. dx<p 
dp = . 

TT 



(C2) 



In the paper, we consider temperatures very low com- 
pared to the Peierls transition temperature, and we take 
1^1 = 1, = 1. This yields Eq. 



APPENDIX D: DERIVATION OF THE 
HAMILTONIAN OF A THREE-DIMENSIONAL 
CHARGE DENSITY WAVE 

In this appendix, we provide a derivation of the Hamil- 
tonian of a three dimensional CDW starting from the 
original Fukuyama one dimensional description. The La- 
grangian of a CDW in a single chain is given byi^ 



C = 



hvp 

47r 



dx 



(Dl) 



Obviously, this Lagrangian describes phason waves prop- 
agating with the velocity u^. vp is the Fermi velocity of 
the electrons forming the CDW and 4> is the phase of the 
CDW. We define an effective mass m* by: 



m 



(D2) 



where Wg is the electron mass. 

In a three-dimensional CDW with screened Coulomb 
interactions, the chains are coupled by a backscattering 
interaction. The resulting Lagrangian reads: 



hvp X - 



dx 



{x,n) (D3) 



-f- ^ J(n, n') / c?a;cos[(/)(a;, n) — (/)(x, n')], 

(n.n') 

where J(n, n') is short-ranged and is given by Eq. (|B6 
Expanding the cosines, 



= 1 
= 1 

and defining. 



cos[0(x, n) — (f){x, n -I- Gy)] = 
(0(x,n) - (/)(a;, n + Gy)) 



V I Pi J,^2 



hvp V. 



47r vl 



J. 



(D4) 



(D5) 



10 



the Lagrangian in Eq. ljD3|) can be rewritten asS^: 



47r 



(D6) 



The sum over lattice sites in the transverse direction can 



C = 



an integral, by writing: 






4^44 J 


[ -I 







(D7) 



The phason dispersion is now Lo{qY = v'^q'^ + VyQy+Vzqt, 
The momentum conjugate to is obtained by the usual 
relation: 



hvf 



(D8) 



yielding the Hamiltonian 



H =. 



Airt 



2 

y_ 

2 ^"V 



(D9) 



From the Hamiltonian Eq. (jDQj, it is straightforward 
to obtain the Debye- Waller factor associated with the 
zero point fluctuations of the phase (f). In the isotropic 
case, Vy = Vz — f0, one finds that (cos ^(a::))T=o ~ 
g-c(me/m where C ^ tt/A is a dimensionless con- 
stant of order 1. Due to the smallness of the ratio 



me/m* ^ 10~^, the zero point motion can be neglected, 
and the kinetic term oc 11^ in Eq. HD9() can be dropped. 
This leads to the Hamiltonian H16|l . 



APPENDIX E: ESTIMATION OF THE 
DISORDER STRENGTH 



Here, we give an estimation of the disorder strength 
D in doped KMoi-xVxOa. We assume a binomial dis- 
tribution of Vanadium impurities on the Molybdenum 
sites. The Vanadium impurities carry an extra electron 
compared to the Molybdenum ions. The resulting charge 
density fluctuation reads: 



5p{v) = ^(a; - ai^a)S{v - Ri,a), 



(El) 



where i is the index of the cell and a is the index of 
the Molybdenum site in a given cell, ai^a = if the 
site is occupied by a Molybdenum ion, and Ui^a = 1 if 
it is occupied by a Vanadium impurity. By construction, 
the expectation value of (5p(r) is zero. We estimate the 
second moment of 5p{v) as: 



5p{v)5p{v') = 

= xil -x)J2 '^(i" - ^^,c.)Sir' - r). 



(E2) 



where we have used the property that 
(x - tTj,Q,)(a; - CTj^/s) = SijSa^f3ix - aj^fi)^. The ex- 
pectation value of ^ i5(r — Ri^a) is simply the number 
of Molybdenum ions per unit volume, leading to the 
formula Eq. l(5T|l . 
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